BROWNIAN MOTION WITH VARIABLE DRIFT 
CAN BE SPACE-FILLING 



TONCI ANTUNOVIC, YUVAL PERES AND BRIGITTA VERMESI 

Abstract. For d > 2 let i? be standard d-dimensional Brownian motion. For any 
a < 1/d we construct an a-H61der continuous function /: [0, 1] — ?> R'' so that the range 
oi B ~ f covers an open set. This strengthens a resuh of Graversen (1982) and answers 
a question of Le GaU (1988). 
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i 1. Introduction 

For d > 2, let Bt be a c?- dimensional standard Brownian motion and / : [0, 1] — )■ Mf^ 
continuous function. We say / is polar for rf-dimensional Brownian motion if, for all x, 

: ^ABt = fit) for some t > 0} = 0. 

+^ . It is well known that Bf does not hit points (see Corollary 2.24 in [B]) almost surely, 

hence constant functions are polar. The Cameron-Martin theorem (see Theorem 1.38 in 
[S] or Theorem 2.2 in Chapter 8 in [B]) tells us that for all functions / in the Dirichlet 
space D[0, 1] (integrals of functions in L^[0, 1]), the path Bt — f{t) has the same almost 
sure properties as the Brownian motion path. Hence functions in D[0, 1] are polar. For 
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CN ■ dimension two Graversen proved the following. 
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2 ■ Theorem A (Graversen [3]). For a// < 7 < 1/2, there exists a 'j -Holder continuous 
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function f : — t- which is non-polar for planar Brownian motion. 



Graversen also conjectured that all 1/2-Holder continuous functions are polar and gave 
partial answer to this conjecture. In 
a similar result for dimensions d > 3. 



' a partial answer to this conjecture. In 15J, Le Gall proved Graversen's conjecture and gave 



Theorem B (Le Gall [5J). Let f : M.'^ be a continuous function. Then f is polar 

for d- dimensional Brownian motion if it satisfies one of the following conditions: 

(i) d = 2 and for all K > there exists a 6 > such that for all s, t with < s < t < 
K and < t - s < S, 
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(ii) d> 3 and for each K > 0, 

,„p I/W-/WL0. (1,2) 

t~s<5 

At the end of [5] Le Gall asked whether for each 7 < l/ii there exist 7- Holder con- 
tinuous functions which are non-polar for (i-dimensional Brownian motion. By a simple 
application of Fubini's theorem, a function / is non-polar if and only if the expected vol- 
ume of the range of B — / is positive. In the following theorem we give a positive answer 
to Le Gall's question and also strengthen Graversen's Theorem. 

Theorem 1.1. For d > 2, let Bf be a d- dimensional Brownian motion. Then for any 
7 < 1/d, there exists a 'j-Holder continuous function f : [0, 1] — ?■ M.'^ for which the range 
of Bt — fit), for t G [0, 1], covers an open set almost surely. 

Note that, in Theorem [Bl the condition for polarity in higher dimensions is stronger 
than for 2 dimensions. We prove the following theorem which shows that a stronger 
condition is indeed necessary. 

Theorem 1.2. For d > 3, there exists a l/d-Holder continuous function / : [0, 1] — )■ M'^ 
that is non-polar for d-dimensional Brownian motion. 



In the proof of Theorem 11.21 we take / to be the standard d-dimensional Hilbert curve. 

2. Hilbert curves 



To prove Theorem 11.11 we will construct space filling curves that remain space filling 
even after being perturbed by Brownian motion. These curves are modifications of the 
standard Hilbert curve construction which is briefly described below, for dimension two. 

2.1. The standard Hilbert curve. Start with the unit square [0, 1]^. In the first iter- 
ation, subdivide the unit square into four sub-squares of side length 1/2 and join their 
centers via three line segments. Choose a direction for the resulting path. In iteration 
n + l, take the path from iteration n, rescaled by 1/2, and put one copy in each of the four 
sub-squares, rotated clockwise by 7r/2, 0, 0, 3n/2 in the bottom- left, top-left, top-right 
and bottom-right sub-square, respectively. Then connect the four copies by adding three 
connection segments. The limiting curve obtained by following this iterative construction 
is called the Hilbert curve. Figure [T] illustrates the first four steps in the standard Hilbert 
curve construction, where we always traverse the path starting at the left bottom square 
and ending in the right bottom square. 

An equivalent description of the Hilbert curve is as follows. For n > 1, at iteration 
n, the unit square is partitioned into a collection of 4*^ squares of side length 2~". Each 
of the 4" squares are labeled by some string ai, . . . , a„ with entries in {0, 1, 2, 3}, which 
determines the order in which the squares are traversed. The main idea is to find an 
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(a) n=l (b) n=2 (c) n=3 (d) n=4 

Figure 1. The first 4 steps in the Hilbert curve construction 

algorithm that recursively determines the order in which sub-squares are traversed given 
the order of their "ancestor" squares. Not picking the right order causes the limiting 
function not to be continuous. The requirement is that for any sequence ai . . . a„ of 
elements in {0, 1, 2, 3} with a„ G {0, 1, 2}, the square labeled by ai . . . a„33 ... 3 is very 
close to the square ai . . . (a„ + 1)00 ... 0. 

Definition 2.1. The traversal order will be described by functions h satisfying: 

(i) h : {0, 1, 2, 3} ^ {(-1, -1), (1, -1), (1, 1), (1, -1)} is one-to-one, 

(ii) M{0,2}) = {(-1,-1),(1,1)}, 
R M{1,3}) = {(1,-1),(-1,1)}. 

Suppose the string ai,...,an denotes the labeling of a square in The traversal 
order of its four sub-squares is given by the function fai,...,a„, which we now construct 
recursively. For the empty string define /0 as 



(/0(O), /0(1), /0(2), /0(3)) = ((-1, -1), (-1, 1), (1, 1), (1, -1)). 



(2.1) 



Algorithm 1 (Standard Hilbert). Given /ai,...,a„_i and a„ define /ai,...,a„ to satisfy con- 
ditions in Definition 12.11 and 

• if a„ G {1,2}, let fa^,...,a„ = /ai,...,a„_i- 

• if a„ = 0, let /ai,...,a„(0) = fa,,...,a„-M, /ai ,...,a„ (l) = /ai ,...,a„_i (3) , 

• if a„ = 3, let /ai,...,a„(0) = /ai,...,a„_i (2), /ai,...,a„(l) = /ai ,...,a„_i ( 1) • 

Now it is not hard to give a more precise definition of the standard Hilbert curve. 
Let O.Oi . . .an ■ ■ ■ denote the 4-ary expansion of t G [0, 1]. We construct functions Hn '■ 

[0, 1] [0, If by setting Hi{t) = ( ^, i ) + ^/0(ai) and for n > 2, 



r 2 

Hn{t) = Hn-i{t) 



n+l 



ai,...,a„_i \'-^nJ 



(2.2) 



Let H{t) = lim Hn{t) be the resulting Hilbert curve. H is well defined, space filling and 

n— >oo 

1/2-Holder continuous (see Chapter 2 in [9]). Hence, by Theorem [Bl H is polar for planar 
Brownian motion. Also note that the Hilbert curve is self-similar, through the mappings 
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from the construction, and measure preserving with Lebesgue measure on both the unit 
interval and the unit square. 

2.2. Generalized Hilbert curves. We now modify this construction in order to obtain 
curves that are space filling and Holder continuous with the exponent slightly less than 
1/2, but which remain space filling after perturbation by planar Brownian motion. The 
curve will be constructed on overlapping squares, using the same algorithm as above, but 
a different parametrization. 

2.2.1. Main idea. Fix a > 1/2. At level 1, we subdivide the unit square into four sub- 
squares of side length 1/2. Cover each of these sub-squares with a slightly larger (centered) 
square of side length a, as in Figure [21 Iterate the procedure. Thus, in step n, each of 
the 4" squares of side length a" are covered by four overlapping squares of side length 
want to find a labeling that determines the order in which these overlapping 
squares are traversed, so that we obtain a continuous function in the limit. Naively, let 
us first try the labeling and parametrization used for the standard Hilbert curve, and to 
time k/A"^ assign one of the 4" squares of side length a" as in the Hilbert algorithm. The 
curve in the first iteration is the same as the standard Hilbert curve and the next two 
iterations are shown in Figure [2J 
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(a) n=l (b) n=2 (c) n=3 

Figure 2. First 3 steps in the new construction with a = .9. The blue 
square is the unit square, red lines are connections between copies of lower 
level iterations, and the gray squares denote typical covering squares at each 
level. 

The limiting function is space filling, but it is not continuous. To see this, first note 
that following the same algorithm as in the Hilbert curve construction leads to eventually 
starting and ending the curve outside of the unit square. In figure O when a = .9, this 
happens already in the third step. For a given a, there will be an e and an integer m 
such that the starting and the ending point of the curve after m*'^ iteration is at least e 
distance, in each coordinate direction, outside the unit square. Following the standard 
Hilbert algorithm the endpoints of the part of the curve in the top-left sub-square should 
converge to the image of 1/2. The same is true for the starting points of the part of 
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curve in the top-right sub-square. However after m*^ iteration these starting points and 
endpoints will be (in each coordinate direction) at distance at least ae from the center of 
the unit square. By symmetry their distance will be at least 2ae. Therefore their limits 
will not be equal which implies the discontinuity at 1/2. 

The main idea in our construction is the following: in each iteration, we assign three 
intervals of time reserved for connecting the four descendents. Choosing the right way 
to assign these connection times, along with the appropriate expansion factor for the 
squares, will allow us to obtain a Holder continuous curve of any exponent less than 1/2. 
The details for this construction can be found below. 



2.2.2. Scaled squares. Let < a < 1 be given and set (5 = which might be negative 
if a< 1/2. 

Analogous to the standard Hilbert curve construction, the set of squares at iteration 
n + 1, denoted by J^n+i, is constructed from J^n by taking J^q = {[0, 1]^} and if [a, a + 6]x 
[b, b + 6] is in J^n, its descendents in J-^+i are 

[a + (i/2 - P)5, a+((i + l)/2 + m x [b + (j/2 - P)5, b + ((j + l)/2 + m, 

for i,j G {0, 1}. We define the partial ordering on U„>oJ>t by setting Qi < Q2 if Q2 is a 
descendent of Qi and extending it by transitivity. Note that the side length of squares 
in J-'n is equal to a"" and if {x, y) is the center of a square in J^n then the centers of its 
four descendents are given by (x.y) + ^a"'{i,j), where i,j e { — 1, 1}. If Qi G J-n then 
^Qi<Q2Q'2 is equal to Qi if a < 1/2 and otherwise is a square with the same center as Qi, 
but with side length (1 + 3^)0;" — 2(i-a) ■ denote the square LiQ^<Q2Q2 by scale((5i). 
In particular if a > 1/2, the union of all squares, over all generations, gives the square 

L l-a' ^ ~ 1-ai ■ 



2.2.3. Generalized Hilbert curve - construction. Let p < 1/4 and p < o; < 1 be given. 
First we construct a Cantor set Cp as follows. At stage 1, divide the unit length interval 
into 7 subintervals out of which we keep four closed intervals of length p. The collection 
of the intervals still kept at stage 1 is 



1-p l + 2p 



2-2p 2 + p 



The removed open intervals will be called connection intervals. At stage n, each remaining 



closed interval of length from C^" is divided into seven subintervals out of which 
four closed intervals are being kept. For example, if [a, a + G C^"^, then we have 

four descendents of this interval 



, a + p" + -(l-p)p' 



n-l 



:A; = 0,l,2,3^cC". 
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Let Cp = n^^Cp be the resulting Cantor set. For every t G Cp, there exists a unique 
sequence {a„}^^ with a„ G {0,1,2,3}, such that t = ^Y^Yl^=i'^np^- Let us call this 
sequence the p- expansion of t. 

Recursively define the sequence of functions G„ : Cp — )■ as follows. Abusing notation, 
denote the p-expansion of t by 0.aia2 . . . and let Gi(t) = (|, |) + ^/0(ai). For n > 2, let 

^^(t) = + ^/a„...,a„_i(a„), (2.3) 

where the functions fai,...,a„ are defined according to Algorithm [H Let G{t) = lim G„(t) 
for all t G Cp and linearly interpolate between these values to obtain a continuous curve 
on [0, 1]. It is an easy exercise to find a formula for G{t) when t ^ Cp, but we will not need 
it here. The function is clearly well defined and continuous. We will show next that G is 
Holder continuous and moreover space-filling (its range covers an open set) for a > 1/2. 

2.2.4. Properties of generalized Hilbert curves. The range of G„ will be exactly the set 
of centers of the squares in For any t G Cp whose p-expansion starts with O.ai . . . a„ 
denote by Dn{t) the set of numbers in Cp whose p-expansion starts with the same n digits. 
If Gn{t) is the center of Q G J-"™ then numbers r G Dn{t) will be mapped by Gn+i into the 
centers of descendents of Q. Therefore Gn+kij) G scale (Q) for all /c > and r G Dn(t), 
in particular after taking limits G{r) G scaAe{Q). 

We claim G is space filling if a > 1/2, more precisely, the image of G covers the entire 
square scale([0, 1]^) = [— 1 + tt^]^- Let U be an open set that intersects the interior 
of the square [—jz^, 1 + jz^]'^- Then U contains a square scale{Q) for some Q G J^n, and 
for a sufficiently large n. As argued before, the image of G will intersect scale(Q) and 
thus U. Therefore, the image of G is dense in [—jz^, 1 -l- jz^]'^ and by compactness, it is 
equal to the whole square [—jz^, 1 -l- jz^]'^- 

We will now compute the Holder exponent for G. First it is important to note that 
if the interval (a, a + i^p") is a connection interval removed at level {n + 1) in the 
construction of Cp then \G{a) - G{a + i^p")| = Ca", where C = ^ if a > 1/2. Also 
note that the endpoints of connection intervals belong to Cp. 

Let t, s in [0, 1], with \t-s\ < ^p^'^ be given. There few cases to be considered. 

(i) If t, s G Cp, then they have common first n digits in their p-expansions. Hence 
Gn{t) = Gn{s) and by the discussion above, there exists a constant Ci, depending 
only on a, such that 

\G{t)-G{s)\ <Cia". 

(ii) If t, s ^ Cp and they belong to the same connection interval (a, a -|- ^—^p^) for 
some k, then 
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for some constant C2 which depends on a and p. Here we have used the fact that 
a/ p > 1 . The same argument holds for s inside a connection interval and t being 
one of the endpoints of this interval. 

(iii) If t G Cp but s ^ Cp, and t < s, let a be the left endpoint of the connection interval 
containing s. Then 

\G{t) - G{s) I < \G{t) - G{a) \ + \G{a) - G{s) | < Ci + C2 a", 

from cases (i) and (ii). Similarly, we use the right endpoint of the connection 
interval if t > s. 

(iv) If t, s ^ Cp but they belong to disjoint connection intervals, and t < s, we let a be 
the right endpoint of the connection interval containing t. Then using cases (ii) 
and (iii), 

\G{t)-G{s)\ < \G{t)-G{a)\ + \G{a) - G{s)\ < Cia" + 2C2a". 
From here we see that G is l^s^. Holder continuous. Note that p can be made arbitrarily 

log p r J 

close to 1/4 and if we let a J, 1/2 and p f simultaneously, we can get a Holder 
exponent arbitrarily close to 1/2. 

Remark 2.2. Choosing different expansion scales for different iterations, one can construct 
a single function that is 7- Holder for all 7 < 1/2. In the n-th level of the construction we 
need to take a„ as the expansion factor for the squares and pn as scaling factor for the 
Cantor set so that | 1/2 and p„ t 1/4- 

3. Proof of Theorem 11.11 {d = 2) 

Let Bf be a planar Brownian motion. Almost surely, there is a (random) e such that, 
for any s < e and any < t < 1 — s, we have |-Bt+s — Bt\ < Sy^s log(l/s) (see Theorem 
1.14 in [6J). The result now follows from the theorem below. 

Theorem 3.1. Let G be the generalized Hilbert curve constructed with a > 1/2. Let 
h: [0, 1] — 7- 6e a function for which there is a constant G such that, for s small enough 
and < t < 1 — s, we have \h{t + s) — h(t)\ < G^ys log(l/s). Then the range of the 
function G + h covers an open set. 

Proof. Pick no large enough such that the assumption on h holds for all s < p""-! g^^^ such 
that /So;" > G^y p^n log(l/p) for all n > uq. Let Qn{p) denote a square of side length a" 
centered at p, and recall the definition of Dn{t) from subsection 12.2.41 By construction, for 
any t G Cp, the square Qn{Gn(t) + h(t)) is covered by the squares Qn+iiGn+i{r) + h{t)), for 
r G Dn{t), even when each of them is shifted by no more than /Ja""*"^. Now by assumption 
on h, the squares Qn+i{Gn+i{r) + h{r)), r G Dn{x) cover the square QniGnit) + h{t)) 
whenever n > uq. Therefore, for all m > n > uq the square Qn{Gn{t) + h{t)) is covered 
by the squares Qm{Gm{r) + h{r)), where r ranges over the set Dn{t). 



8 



Let U be an open set that intersects Qno{Gno{t) + h(t)) for some t E Cp. By the 
discussion above, one can find m large enough and r G Dngit) such that the square Q 
with the center at Gmij) + h{r) and side length equal to C \J p'^m\og{l / p) + 2{i^a) ^^^^ 
?7. If s G Dfn{r), then G'(s) + lies in Q and thus in U. This implies that the image 
of G + /i is dense in the square Qno{Gno(t) + h(t)) and, by compactness, it covers it. □ 



Remark 3.2. As mentioned in Remark l2.2[ by taking different a„ and p„ at different stages 
in the construction of the Hilbert curve, we can obtain a generalized Hilbert curve which 
is 7- Holder continuous for all 7 < 1/2. Actually such curves can be constructed in a way 
so that they satisfy the two-dimensional case of Theorem 11.11 To do this we only need to 
ensure that Brownian motion does not destroy the overlap of squares at different levels 
and the condition for this is that for large enough n we have 



i=l 



\ 



n log ( jj , , 

1=1 1=1 



where Pn+i = ^""^^ Such sequences q;„ and p„ can be constructed with q;„ I 1/2 and 
Pn^ 1/4:. For example take a„ = and p„ = ^ 

4. Construction of a reverse Holder continuous drift 

In this section we describe an alternate construction of a function which added to 
standard Brownian motion will give an almost surely space-filling curve. While its con- 
struction is slightly more complicated, this function has the advantage that it is reverse 
Holder continuous. 

Definition 4.1. We say that a continuous function /: M — > M'^ is reverse a-Holder 
continuous, for < a < 1, if for some C > and any open interval / of length |J| we can 
find s,tel such that - f{t)\ > 

The notion of reverse Holder continuity function is closely related to the geometric 
properties of its graph. For example graphs of functions which are both a-Holder and 
reverse a- Holder continuous (like the function we will construct in this section) have Haus- 
dorff dimension bigger than 1, see Theorem 4 in 0. Clearly, the previously constructed 
generalized Hilbert curve is not reverse Holder continuous since it is linear outside the 
Cantor set C^. 

Let a > 1/2 be given and recall that /3 = {2a — l)/4. Let S* be a set of positive integers 
such that Ylies'^^ ~ i^- This is the set of exceptional times when our construction 
deviates from the standard Hilbert algorithm. One way to construct this subset is to 
define Ai := and, at step i, include in S the smallest integer j, larger than all integers 
already in S, which has the property < Aj and set A^+i = Aj — . We call the string 
ai, . . . , a„ 0-reverse {3-reverse) if ai, . . . , a„ ends in exactly k O's {k 3's) for some k with 
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k + 1 & S. This means an-k+i = ■ ■ ■ = = 7^ o-n-fc, and similarly for 3's. For a string 
that ends in several consecutive O's (3's) we define its prestring to be the string obtained 
by deleting the maximal block of consecutive O's (3's) at the end. 

Algorithm 2 (Alternate generalized Hilbert). Define fm by (12. ip . Let be given 

and let ai, . . . ,0^ be its prestring. Given fai,...,a„^i and a„ define fai,...,a„ to satisfy the 
conditions in Definition 12.11 and 

• if a„ e {1,2}, let /ai,...a„ = /ai,...,a„_i- 

• if a„ = 0, let /ai,...,a„(l) = /ai,...,a„_i(3) and 

- if ai, . . . , a„ is 0-reverse, let /ai,...,a„(0) = /ai,...,a„(2). 

- if ai, . . . , On is not 0-reverse, let /ai,...,a„(0) = /ai,...,a„(0). 

• if a„ = 3, let /ai,...,a„(0) = /ai,...,a„_i (2) and 

- if ai, . . . , a„ is 3-reverse, let /ai,...,a„(l) = /ai,...,a™(3). 

- if ai, . . . , a„ is not 3-reverse, let /«!,...,«„ (1) = /ai,...,a^(l)- 

Forn > 1 let Q„ = {^4"" : A; G Z, < < 4"}. For every < t < 1 consider its 4-ary 
expansion t = 0.aia2 .... Now define Go to be the constant function mapping into the 
center of the unit square [0, 1]^ and for n > 1, let the sequence of functions Gn be given 
by recursion (12. 3 p . This is a Cauchy sequence of functions and so finally define G as their 
limit. 



(a) n=l (b) n=2 (c) n=3 

Figure 3. First 3 steps in the new construction with a = .9. Note that 
the curves start and end inside the unit square, drawn in blue. 

Space filling property of this function is obvious. 

Now we prove that the constructed function is both '"i^g 4"'' -Holder and reverse ^"f^g 4°^ - 
Holder continuous. For the reverse Holder continuity notice that any interval / C [0, 1] 
contains a 4-ary interval J, of the form [k4:~"', (/c+l)4~"], of length at least |/|/8. Moreover 
any 4-ary interval J of length 4"*^ is mapped by G onto a square of side length 2(i-a) ~ 
_L_|j|iog(i/a)/iog4_ pipi^ s,t e J so that G{s) and G{t) are vertices of one side of the 
square G{J). We see that 

\n(Q\ — niiW — I T|log{l/a)/log4 ^ " I r|log(l/a)/log4 

^Wl-2(i_«)l^l -2(l-a)'^' 



10 

If s G Qn and t is any number whose 4-ary expansion agrees in the first n digits with s, 
then \G{t) - G{s)\ < ^(f-a) - Take any < ti < < 1, and n such that t2 - < and 
gi, g2 e Q„ such that <ti-qi < 4"", for i = 1, 2. Then we have \G{ti) - G{qi) \ < ^(fza) 
forz = 1,2 and \G{t2)-G{ti)\ < ^ + \G{qi)-G{q2)\. Since either gi = gs or ga = gi+4-" 
we see that it suffices to prove the Holder continuity for points in IJnQ"- 

Let s, t G Q„ with t = 5 + 4"". If the 4-ary expansion of s does not end in 3 then s and t 
have equal first n — 1 digits and so \G(t) —G{s) \ < -^^^^"^ done. Otherwise let 

the 4-ary expansion of s be given by s = O.ai . . . 0^3 ... 3 with k = n — m threes in the end 
and am 7^ 3. Then obviously t can be written in base 4 as t = O.ai . . . 0^-1 (ctm + 1)0 . . . 0. 
Define si = O.ai . . . 0^3 ... 3 with / threes in the end, so in particular Sk = s. Without 
loss of generality assume /ai,...,a™_i = h- Denote the point Gm-i(t) = Gm-i{s) by p. Let 
us verify the claim for = 1 since cases G {0, 2} are similar. 

If am = 1, then fa^,...,a^ = h = /ai,...,(a™+i) and, relative to p, the coordinates of Gm{,s) 
and Gm(t) are (— °™ , ""^ ^ ) and ( ""^ ; ""^ ^ ) ; respectively. Now, for any I > 0, the 
point Gm+iisi) is constructed from Gm+i-iisi-i) by shifting it by in the lower 

right direction, unless / G S when we have to shift the point by the same distance in the 
upper left direction. From the definition of S it is clear that Gm+i{si) will converge to p 
as / — > 00 and moreover that \Gm+iisi) — p\ < V^a™" X]j>/ — i^'^"'^'- particular 
for / = k this implies that 

\G{s)-p\ < \G{s)-G4s)\ + \Gn{s)-p\ < ^^^"'^ = l±M(t_,)iog(i/«)/iog4_ 

Similarly, we argue that Gm+iif) also converges to p which implies G{t) = p and plugging 
this into (14. ip gives the wanted bound. 

Adding G to planar Brownian motion results in a space filling process. The proof of this 
fact is analogous to the proof of Theorem 11.11 when d = 2, hence we do not include it here. 
In d-dimensions {d > 3), one can also construct a continuous curve with Holder exponent 
arbitrarily close to l/d, with the same properties as G, by changing the traversal algorithm 
given in the d-dimensional Hilbert curve construction. Such a higher dimensional curve, 
analogous to G, can be described using a modification of Butz's algorithm from [2]. 

5. Proof of Theorem 11.11 {d > 2) 

To prove Theorem 11.11 in dimensions higher than 2, we construct generalized Hilbert 
curves in the same way as in two dimensions. There are algorithms that generate stan- 
dard d- dimensional Hilbert curves in dimensions higher than two, which can be realized 
through the labeling algorithm we described for two dimensions. One example is the 
Butz algorithm, see [T] and [2] . These algorithms are completely determined by functions 
fai,...,a„ where G {0, 1, . . . , 2'^ — 1}. Throughout this section /ai,...,a„ will stand for these 
functions. The Hilbert curve is defined as the limit of the functions defined recursively 
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by formulae analogous to (12. 2p . As in two dimensional case, the Hilbert curve is again 
measure preserving. We do not put here the explicit algorithms for generating functions 
/ai,...,a„ as they are relatively complicated. See also [9j for a geometric description of space 
filling curves in higher dimensions. 

Let p < 2"*^ be given and define Qp^d by iteratively removing, at step n + 1, exactly 
2'n.d(^2'^ _ (open) connection intervals of length ^^^rr^p" and keeping 2^""^^^'^ intervals of 
length p"+^ As in the 2-dimensional case, every t G Qp^d has a unique representation 

oo 

t = K^a„p", a„G{0,l...,2'^-l}, 

n=l 

with K = (-2dli)p - We will call 0.aia2 . . . the (p, d) — expansion of t. Fix an a > 1/2 and 
define functions Gq to be the constant function mapping to the center of the unit cube 
[0,1]'^ and recursively by f l2.3p . with the understanding that the traversal algorithm 
and hence the functions fai,...,an correspond to the d-dimensional Hilbert curve. Then 
define as the limit of functions G^. Again, for t ^ Cp d, the values G'^(t) are defined 
by interpolation. The resulting curve is clearly well-defined and continuous. Holder 
continuity and space covering property for follows by the same argument as in two 
dimensional case. The Holder exponent is again We leave the details to the reader. 
Taking p close to 1/2'^, and a close to 1/2, leads to continuous functions with Holder 
exponents arbitrarily close to 1/d. 

Proof of Theorem \l.l\ (d > 2). As in the two dimensional case, there exists a constant 
C > such that, almost surely, for every small enough s > and all < t < 1 — s, 

\Bt+s - Bt\ < C^/slog{l/s). (5.1) 

Suppose no is large enough such that condition (15.11) holds for s < g^^^ g^j.]^ that, 

for all n > no, 

/3a" > C^p^nlog{l/p). (5.2) 

For t G Cp^d with (p, (i)-expansion starting with O.ai . . . a„ define Dn{t) as the set of 
numbers r G Cp^^ whose (p, (i)-expansion starts with the same digits. Take the box Q^q of 
side length a"" centered at Bt — Gf^^it). Just as in the proof of Theorem 13. H it is easy to 
see that the 2^ boxes of side length a"'°~^^ centered at Bj. — G'^^^_^_-j^{r) , for r G Dn^it), cover 
Qno- Iterating, we see that for all m > 1, the box Qn^ is covered by the 2""^ boxes of side 
length a"""*"™ centered at Br — G^y_,_,^(r), for r G DnQ^mif)- Then given an arbitrary point 
X G Qno; one can find a sequence of boxes {Qno+m}m>i, of side length centered as 

above, that contain x. From here it easily follows that x is in the image oi B — G"^, which 
completes the proof. □ 

Remark 5.1. Let G be a generalized ci-dimensional Hilbert curve with a < 1/2. These 
cases are also interesting since B — G is not space fiUing and moreover, for p < a^, we 
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can compute the Hausdorff dimension of its range to be exactly max (^dy^^j^,2^. To 
prove this notice that since [0, l]\Cp^rf is a countable union of disjoint intervals, Hausdorff 
dimension of {B — G){[0, 1]\C p^d) is equal to 2. Therefore it is enough to prove that 
{B — G){Cp^d) has Hausdorff dimension equal to d^^^j^- Denote with Qn{p) the cube 
of side length a", centered at p. It is easy to check that there is a 7 > such that, for 
all t G Cp^d, the cubes Qn+i{Bt — Gn+i{r)) for r G Dnit) will be disjoint and contained 
in the cube Qn{Bt — Gn{t)), even if each of the smaller cubes is shifted by no more than 
7a"'"'"^. Now knowing the modulus of continuity of Brownian motion and using the same 
argument as in the proof of Theorem ll.il (and the fact that p < a^) we see that for some 
(random) uq and all n > uq cubes Qn+i{Br — Gn+i{r)) for r G Dn(t) will be disjoint 
and contained in the cube Qn{Bt — Gn{t)), even after being shifted by no more than 
7Q;'^+^/2. As a consequence, for any m > hq and r G Dm{t) every two of the points 
Bj. — Gmir) are at the distance greater than 7q;™/2. This immediately implies that the 
Hausdorff dimension of the image {B — G){Cp^d) is almost surely bounded by d ^J^^J^^ . 
For the lower bound we can use the Mass Distribution Principle (see 4.2, Chapter 4 in 
[3]). For a fixed t G Cp^d we define a measure on the part of {B — G){Cp^d) which is 
contained in Qno{Bt — Gnoit)). The measure is simply defined by assigning mass 2"'^'^ 
to each of the disjoint cubes Qno+k{Bs — Gna+k{s)) for s G Dn^it)- Then any ball of 
radius 7a'^°"'"'"/2 is assigned no more than 2"'"'^ mass and the claim follows by the Mass 
Distribution Principle. 

Therefore, for any 2 < \ < d and any /t < 1/A, we can find a generalized d- dimensional 
Hilbert curve which is K-Holder continuous and such that almost surely the range of 
Brownian motion with this drift has Hausdorff dimension equal to A, by taking a = p'^ = 
2-d./>' in our construction. 



6. Proof of Theorem 11.21 

We will prove Theorem 11.21 by showing that for > 2 the range of dimensional 
Brownian motion with drift given by the standard rf-dimensional Hilbert curve hits points 
with positive probability. The key ingredient in the proof is the following lemma which 
shows that, for any a > 1/d, the Hilbert curve satisfies the reverse Holder inequality for 
all points from a set of large measure. 

Lemma 6.1. For d > 3, let H : [0, 1] — ?■ [0, 1]^ denote the standard d- dimensional Hilbert 
curve. Then for any a > 1/d and any e > 0, there is a set S C [0, 1] of Lebesgue measure 
greater than 1 — e and a G > such that \H{s) — H{t)\ > G\s — t|" for all s,t & S. 

Proof. As usual is the set of dyadic cubes used in the construction of the Hilbert curve. 
Define Wn := UQeJ"„ '^Q- ^'^^ 6 > we define B{Wn, 6) as the set of points in [0, 1]°' at 
distance no more than 6 from Wn- 
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First we show that the volume of the set B{Wn,2~"~"^) is bounded from above by 
2d2~'^. We prove this by induction. For n = 0, this is trivial since Wq is just the 
boundary of the original cube whose {d — l)-dimensional volume is equal to 2d. When 
m = 0, the claim is also trivial since B{Wn,2~"'~"^) is equal to the whole original cube. 
In the general case, notice that the intersection of B{Wn, 2~"~™) with a cube from J-'i 
is equal to -B(l^„_i, 2^"^™+^), scaled by a factor of 1/2 and translated accordingly. The 
induction assumption implies that the volume of this intersection is bounded from above 
by 2d2-'^2-™. Summing over all cubes in J-'i proves the claim. 

For c> define K := U„>o B{Wn, 02""'='") and set S := H-\K^). Since the volume of 
K is less than cJ2n>o'^~''^'^'^~^^ and H is measure preserving, we can choose the constant 
c to make the Lebesgue measure of S arbitrarily close to 1. 

To check that S satisfies the desired property, take s,t E S and let n be the largest 
integer so that H{s) and H{t) are contained in the same box in Tn-\- Then |s — t| < 
2-d{n-\)^ but \H{s)-H{t)\ > 2c2-"^", since H{s) and H{t) are in and inside different 
boxes in J-'n- This proves the claim with C = 2^~'^c. □ 

Proof of Theorem M.B, Let 1/d < a < 2{d+i) — 1/2 be given and let H denote the regular 
d- dimensional Hilbert curve. Let J be a set of positive measure with \H[s) — H{t)\ > 
C\s — t\°' for all s,t G J, for some constant C > 0, as in Lemma 16. 1[ Without loss of 
generality we can assume J is closed. 

Fix 6 > and x G Mf^, and set Jg = [6, 1] fl J. We let be the amount of time in 
Js that the range of function t ^ Bt — H{t) spends in := B{x,e), the ball of radius 
e around x. We will apply the second moment method to T^. Using Fubini and the fact 
that the density of (i- dimensional Brownian motion is bounded from below on when 
t > 6, we have 

E(T,) = [ f>{Bt- H{t) eB,)dt>he\ (6.1) 

for some constant h depending only on 5 and d. To find the second moment of T^, apply 
Fubini to obtain 

E{T^) = 2 [ F{Bs-H{s)eB,) [ F{Bt - H{t) e B,\B, - H{s) e B,)dtds. (6.2) 

Jjs Jjsr\[s,i] 

By bounding density for s G J<5 we have P(-Bs — H{s) G B^) < CiVde'^, where Ci depends 
on 6 and denotes the volume of the unit ball in d dimensions. Thus we only need to 
estimate the inner integral in (16.21) . Define the set 



Js = {te J -.tys, \H{t) - H{s)\ > 3e}. 
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The inner integral in fl6.2p can be separated into three parts. Let 



Ai = I ¥{Bt- H{t) e B,\Bs- H{s) e B,)dt 

[s,i]n{Js\Js) 



A2 := / - H{t) e B,\Bs - H{s) G B,)dt 

As:= I F{Bt- H{t) E B,\Bs~ H{s) E B,)dt. 

ifs.s+f4e/c)i/QinjAnJ» 



[s,s+(4e/c)i/Q]nJ4nJ. 

Because H is measure preserving, the Lebesgue measure of J\Js can be bounded above 
by Vd3'^e'^. Thus, using a trivial bound for the probability inside the integral, we have 

A, < v,3''e''. 

For the second integral, notice that the distance between the balls B{x + H{s),e) and 
B{x + H{t), e) is at least C(t — s)" — 2e > C{t — s)" /2. Therefore, by bounding the density, 

A2< / Vde'' ; ^ U t<C2Vde'', 

Js+(4e/c)^/- (2vr(t - s))'='/2 

for some constant C2, where we use the fact that limrj^o exp(— rP)r'^ < 00 for p < and 

any q. Finally, if both Bs — H{s) and Bt — H{t) are in B^ and t E Jg, then it must hold 

that \Bt — Bs\ > e and so 

A3 < / P ' * ' > , dt 

l'S+{4e/Cy/°' 2(7' / 4 \ 3/(2a) , 

< / c'^^^dt = ^(^) 

y g 6 o V O / 

where the second inequality follows by Markov inequality. Since a < 2{d+i) ' follows 
that A3 < C3e'^, for some constant C3. 

Now summing all the bounds, we have that K{T^) < Ke^'^ , for some K > and by 
Paley-Zygmund inequality we obtain 

¥(Bt - Hit) = X, for some t) > limPfTe > 0) > hm '\ > — . 

□ 

Theorem 11.21 tells us that even after the Brownian perturbation the standard Hilbert 
curve has positive volume. It would be interesting to see whether it also remains space 
filling. 

We end the paper by giving a proof, different from the one in [5], of Le Gall's result 
that any function satisfying (11.21) is polar. It is based on the fact that images of such 
functions have zero volume. In the following 15*1 will denote the diameter of a set S. If 
/: ]R+ — R"' satisfies (11. 2p then for any e > we can find a 5 > such that for any set 
5 C M of diameter 15*1 < 5 we have < e|5'|^/"'. Covering an interval / of length I by 

intervals Ii, ... In of length l/n for n large enough, we obtain a cover f{Ij), 1 < j < n for 
the image /(/) which satisfies ^ . \f{Ij)\'^ < el. Therefore the d-dimensional Hausdorff 
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measure of /(/) is zero. Since Brownian motion is almost surely a-Holder continuous for 
any a < 1/2 we have that B — f satisfies the same condition (11 ■2p and therefore its image 
has volume zero almost surely. Now the result follows by an application of Fubini and 
standard arguments. 
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